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We present a method for determining correlations in a gas of indirect excitons in a semiconductor 
quantum well structure. The method involves subjecting the excitons to a periodic electrostatic 
potential that causes modulations of the exciton density and photoluminescence (PL). Experimen¬ 
tally measured amplitudes of energy and intensity modulations of exciton PL serve as an input 
to a theoretical estimate of the exciton correlation parameter and temperature. We also present 
a proof-of-principle demonstration of the method for determining the correlation parameter and 
discuss how its accuracy can be improved. 


I. INTRODUCTION 

Indirect excitons (IXs) in coupled quantum well struc¬ 
tures (CQW)^’^ is a model system for exploring diverse 
physical phenomena including pattern formation/ spon¬ 
taneous coherence and condensation/^® transport/^^° 
spin transport and spin textures/^ and localization- 
delocalization transitions.The CQW consists of a 
pair of parallel quantum wells separated by a narrow tun¬ 
neling barrier and the IX is a bound state of an electron 
and a hole confined in the opposite quantum wells. Inter¬ 
actions play a key role in the physics of IX systems. In 
the first approximation, the exciton-exciton interaction 
potential is given by 


v{r) 



-\- (P 


( 1 ) 


where d is the center-to-center separation of the two wells 
in the CQW and k is the dielectric constant of the semi¬ 
conductor. At r ^ d this potential has the form of the 
dipole-dipole repulsion, v{r) ~ e^d?/nr^. A hallmark of 
interacting IXs is the increase of its photoluminescence 
(PL) energy E with density n, which has been known 
since early spectroscopic studies of CQWs.^®’^® At small 
n, where the interactions are negligible, this energy ap¬ 
proaches the creation energy E{0) of a single IX. As n 
increases, a pronounced PL blue shift 


/\E{n) = E{n) - E{0) > 0 (2) 


develops. The physical origin of this energy shift is the 
renormalization of the IX dispersion. The bare dispersion 
is given by 


+ £: ■ <=" 

where k is the exciton momentum, see Fig. 1(a). (In this 
paper the bare quantities are denoted with bars.) The 
renormalized dispersion ek is the solution of the equation 
= £k + X'(ek), where E' is a real part of the self-energy 

( 4 ) 


The simplest theoretical model of the self-energy is the 
Hartree approximation. It assumes that E is equal to the 
energy-independent, real constant 


Tic = nv, V = 


v{r)d^r = 


47re^d 
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( 5 ) 


so that the PL energy shifts by the same amount: 

47rp^ 

AE = Ec = - nd. (6) 

Kj 

Equation (6) is colloquially known as the “capacitor” for- 
mula^^ because it is similar to the expression for the volt¬ 
age on a parallel-plate capacitor with surface charge den¬ 
sity ±en on the plates. The capacitor formula provides 
a qualitative explanation for the observed monotonic in¬ 
crease of AE with photoexcitation power. However, an¬ 
alytical theory beyond the Hartree approximation^®^®® 
and Monte-Carlo calculations®®’®^ suggest that the ca¬ 
pacitor formula significantly overestimates AE{n). The 
same conclusion follows from the analysis of the small-n 
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FIG. 1. (Color online) (a) Bare and renormalized exciton 
dispersions (thin and thick solid lines, respectively; both are 
schematic). The dashed line is the light dispersion; the shaded 
region to the left of it is the radiative zone, (b) The range 
of values for exciton self-energy E (shaded). The upper solid 
line is the Hartree (capacitance) self-energy Ec to which E is 
normalized in this plot. The lower solid line is a lower bound 
estimated from the self-energy of an exciton crystal with pa¬ 
rameters d = 1.2Oe, fle = 10nm, and = 2me suitable for 
GaAs CQW studied in this work. A hydrogenic exciton wave- 
function b(’'eh) = exp(—reh/a) was assumed with a = 2.9 Oe. 


E(e) = E'(e)-hiE"(£). 
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FIG. 2. (Color online) Illustrations of indirect exciton phases: 
(a) crystal (b) liquid with short-range correlations (c) hypo¬ 
thetical uncorrelated liquid or gas. The blue (red) patches 
represent electrons (holes) confined in separate 2D layers. The 
bottom row shows the schematic plots of two-body correlation 
functions G(r) = n~^ (n(r)n(O)) — S(r) for each of the phases. 


and the large-d limits, where IX should form, respec¬ 
tively, a correlated liquid and a crystal, see Figs. 2(a) 
and (b). In such phases the IXs avoid each other, which 
lowers their interaction energy per particle as well as E' 
compared with the uncorrelated gas-like state assumed 
in the Hartree approximation. Fig. 2(c). 

For comparison of various theories with experiment 
and with one another, we define the dimensionless corre¬ 
lation parameter^^ 


where 


7 = vi/S.E/ n , 


= Wx / 27rli^ 


( 7 ) 


( 8 ) 


is the bare exciton density of states (DOS) per spin and 
rux = TOe + TUh [Eq. (3)], We, and m^,, are the ef¬ 
fective masses of, respectively, excitons, electrons, and 
holes. The capacitor formula ( 6 ) predicts the density- 
independent correlation parameter 


7 c = viv = 


2d rUe + mh 


me 


( 9 ) 


which is about 7 c « 7 for the GaAs CQW structures 
studied in Ref. 22 and the present work. Here Ug = 
h? k/{ mee^) is the electron Bohr radius. The interaction 
part of the self-energy of a classical crystal gives a lower 
bound on E' and thus 7 . However, this leaves one with 
a large uncertainty, see Fig. 1(b). A simple remedy for 
the inaccuracy of the Hartree approximation is to replace 
the coefficient v in Eq. (5) with a smaller number t. This 
yields AA = T,' = nt and 


7 = . 


( 10 ) 


The problem remains how to reliably evaluate the phe¬ 
nomenological parameter 7 . The solution of this problem 


can provide both the test of available theories and a con¬ 
venient density calibration tool in experiments. 

In our previous work^^ we outlined a method that al¬ 
lows one to estimate 7 experimentally. The key idea of 
the method is that the same correlation parameter deter¬ 
mines both the PL shift and the ability of IXs to screen 
an external potential perturbation. In other words, the 
screening efficiency of the IX gas is related to 7 . In the 
experimental part of this work, we infer the screening 
properties of the IX gas from variation of its PL energy 
and intensity as a function of coordinates. To do so we 
employ a periodic external potential 


U{x) = ^Ui{l- cosqx) 


( 11 ) 


created with the help of interdigitated gate electrodes. 
The “depth” Ui of this potential is controlled by the gate 
voltage. As the IXs attempt to screen the potential, their 
density becomes periodically modulated (Fig. 3). In our 
previous work,^^ we measured modulation of the IX PL 
intensity and energy in such a device. Analyzing such ex¬ 
periments using a simplified theoretical model, we arrived 
at a rough estimate 7 ~ 2 , which is about one-third of 7 c. 



X ((im) 


FIG. 3. (Color online) Simulated PL energy and intensity of 
indirect excitons in a lattice potential. (a,b) Lattice potential 
[dashed line in (a)], PL energy [solid line in (a)], and exciton 
density (b) for C = 3 meV (measured with respect to eo). (c,d) 
PL intensity (c) and energy (d) under the same conditions for 
NA = 0.4. The depth of the lattice Ui = 4meV, temperature 
T = 4.2 K, correlation parameter 7 = 7 . 
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FIG. 4. (Color online) Localization-delocalization crossovers. 
Indirect excitons are localized in region L, localized along 
X but delocalized along y in region h^Qy, and delocalized 
in both directions in region D. The last of these is suit¬ 
able for measuring 7 . The squares are data from our earlier 
work.^^ The lines are empirical fits AE = (0.8^ -I- 
and Ai? = 0.8 (in units of meV). 


This analysis requires large enough C7; > T ~ 0.8 meV 
to ensure that the lattice potential dominates over the 
random potential of disorder. In addition, sufficient IX 
density producing IS.E higher than Ui and T is required 
in order to overcome localizations effects of both these 
potentials (Fig. 4) as detailed in Sec. II. In the present 
work we expand these earlier efforts into a fully fledged 
method for determining correlations in a gas of IXs from 
the measured modulations of exciton PL in a periodic 
electrostatic potential. Our refined model takes into ac¬ 
count optical resolution effects neglected in Ref. 22. 

It is worthwhile to note that a conceptually similar 
technique has been previously employed for cold atoms 
in magnetic traps. The equilibrium density profile of the 
atomic gas in the trap can be understood as a result of 
screening the trapping potential by the atoms. Further¬ 
more, it has been shown that such density profiles can 
be measured optically and used to determine the equa¬ 
tion of state of the gas.^®’^® By analogy, one can think 
of our external potential U{x) as a periodic array of in¬ 
dividual traps. This geometry offers certain advantages 
compared to single traps. One is the redundancy of the 
experimental data, which reduces stochastic errors of the 
measurements using uniform arrays. Another advantage 
is the simplicity of modeling the optical resolution effects 
in a periodic geometry. 

The remainder of the paper is organized as follows. In 
Sec. II we define the model and derive the main equations 
of the theory. In Sec. Ill we present a proof-of-principle 
demonstration of the method for determining the corre¬ 
lation parameter. In Sec. IV we discuss how to improve 
the accuracy of the method. 


II. THEORETICAL MODEL 

Our objective is to describe PL characteristics of an 
exciton gas subject to a periodic external perturbation, 
Eq. (11), which makes the exciton density modulated. 
However, it is instructive to start with the case Ui = 0 
where the exciton gas is uniform and the in-plane mo¬ 
mentum k is still a good quantum number. 

A. Uniform exciton gas 

The PL intensity spectrum /(w) in our model is given 
by the formula 

/(w) = H/(a; -/r)Ao(w), (12) 

where i? is a constant related to the exciton-photon ma¬ 
trix elements, 

m = (13) 

is the Bose-Einstein distribution function, /3 = 1/T is 
the inverse temperature, and /i is the exciton chemical 
potential. The optical density (OD) Ao{uj) is the k = 0 
part of the spectral function 

A(w) = -2Im<-. (14) 

w - Ek - Xk 

Derivation of Eq. (12) can be found in earlier works 
on the subject.When interactions and disorder are 
present in the system, the functional form of OD may 
be complicated and qualitatively different in the canon¬ 
ical phases of bosonic matter: normal fluid, superfluid, 
and Bose glass.We focus on the intermediate density 
range, 

’T-d > n > y = — , (15) 

t 91 

in which excitons are expected to form the normal fluid. 
Parameter g in Eq. (15) is the number of exciton spin 
flavors (g = 4 in GaAs) and 

nd = gviT (16) 

is the density scale corresponding to the onset of quan¬ 
tum degeneracy. We assume that the interactions and 
disorder shift the exciton energies by the k-independent 
amount = nt, so that 

ek=£k + nt, (17) 

and that the exciton scattering rate |X"| is small com¬ 
pared to the temperature T. Under these simplifying 
assumptions, the OD as a function of uj can be approxi¬ 
mated by a single sharp peak centered at the renormal¬ 
ized exciton band edge Eq = £0 + nt. 
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We restrict our analysis to the two most basic charac¬ 
teristics of the PL spectrum, the total intensity and the 
average energy, which we define using the moments of 
the lineshape: 


=yw'"/(w)^. (18) 

0 

Specifically, the PL intensity is the zeroth moment and 
the energy is the ratio of the first and zeroth moments: 

/ = E = I^^'>/. (19) 


It is easy to find the relations among the key quantities 
in this model: 


n = —rid In 
7= ^ 


1 _ g/5(M-eo) 


II = fj, + nt, 


g^feo-M) _ 1 

E = So = £o + nt . 
Equation (20) follows from 


= B - l) 


( 20 ) 

( 21 ) 

( 22 ) 


OO 

= 9 J - ii)du} 


(23) 


where vi{uj) = vi 6 {lo — £q) is the renormalized DOS. [By 
virtue of Eq. (17), the renormalization is simply a shift 
of the argument by nt.] The first equality in Eq. (21) 
comes from Eq. (12) and the normalization condition 

OO 

J Ao{u;)^ = l (24) 


for the OD. Equation (22) establishes the correspondence 


E{0) = eo , AE{n) = nt (25) 

with the notations introduced in Sec. 1. While the density 
dependence of E is linear, that of fi is not, see Fig. 5(a). 
The nonlinear relation between /i and p, at several dif¬ 
ferent 7 is illustrated by Fig. 5(b). This relation charac¬ 
terizes screening properties of the exciton fluid discussed 
next in Sec. IIB. 


B. Periodically modulated exciton gas 

Let us examine the effect of the external potential U(x) 
[Eq. (11)] on the exciton normal fluid. If U{x) does not 
vary too rapidly, then the semiclassical Thomas-Fermi 
approximation (TEA) is valid. Within the TEA the self¬ 
energy is considered to be position-dependent, S(a:) = 
n{x)t + U{x). It plays the role of the effective potential 





p (meV) 


FIG. 5. (Color online) Exciton density, chemical potentials, 
and effective dielectric constant, (a) Density n as a function 
of the (renormalized) chemical potential /i for (left to right) 
7 = 0.7, 3.5, 7.0. (b) The bare chemical potential as a function 
of the renormalized one for the same 7 (solid lines), (c) Inverse 
e* as a function of fj,. The dashed lines in (b, c) depict the 
7 —>■ 0 limit. T = 4.2 K for all the curves. 


acting on the excitons. The renormalized band edge £q (x) 
is also position-dependent, tracking E(x): 

£ 0 ( 2 ;) = eo + S(a;) = £0 -b n{x)t + U{x). (26) 

The local chemical potential ii{x) is defined such that 
Eq. (20) is unaltered. Finally, the occupation factors in 
Eqs. (12), (23), etc., change from /(w — ^) to /(w — C), 
where 


c = Ii{x) + U{x) (27) 

is the electrochemical potential. One should not confuse 
C, which is cc-independent in equilibrium, with the spa¬ 
tially varying bare p{x) and renormalized /r(a;) chemical 
potentials. By examining the last two we can analyze the 
screening of the external potential by the IXs. Indeed, 
from the above formulas we deduce 

£o(a:) = (£0-f C) -/i(a;), (28) 

lj.{x) = C - U(x). (29) 

Therefore, the plot of £0 vs. U can be obtained from 
the plot oi p vs. /r by inversion and a shift of the axes. 
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This construction is illustrated in Fig. 6 (b) and (e). The 
derivative 


€* 


dU_dfx_ dn 
dYj dp, ^ dp 

1 + (?7 - l ) 


1+37 


exp 


/ £o - gp \ 

V 91T ) 



(30) 


can be regarded as an effective dielectric constant. It 
describes how much variations of the effective potential 
S(a:) are reduced compared to those of the external one, 
U{x). According to Eq. (30), at fixed n stronger inter¬ 
action (higher 7 ) leads to stronger screening (larger e*). 
Conversely, at fixed gp, higher 7 result in lower e*, even¬ 
tually approaching the interaction-independent limit 


e* 


= 1 + 


gp — gp 
T 


(31) 


In order to account for the optical resolution, the convo¬ 
lution with the point-spread function (PSF) of the opti¬ 
cal system must be added. For an ideal circular lens it is 
given by^® PSF(a;) = \Ii{2Qx)/{'kQx^), where 

Q = NA X fcp (36) 

is the largest momentum admitted by the lens, NA < 1 
is the numerical aperture, and 

fcp = w / lie (37) 

is the photon momentum in vacuum. The plot of the 
PSF for NA = 0.24, which is typical of our experimental 
setup, is shown in Fig. 7(a). It implies that the char¬ 
acteristic width of the PSF is comparable to the lattice 
period K = 2tt/ q = 2 /rm, so that accounting for the op¬ 
tical resolution effects is important. Having to evaluate 
the Struve function^® Hi( 2 ;) makes working with the PSF 
inconvenient. Instead, we can write the desired convolu¬ 
tion as the Fourier series. 


If the fixed parameter is 7 , the screening evolves from 
negligible (e, —)■ 1 ) to perfect (e* —)• 00 ) as either n, gp, 
or fi increases. The last case is illustrated by Fig. 5(c). 
Small changes in n and U, which are related by 

1 — 

6n{x) = —6U{x) -, (32) 


Iix,uj) = E Ik{oj) cos kx , 


(38) 


where fc = 0, 3 , 2^, 83 ,... are the Fourier momenta and 

4(w) = (2 - 6k,o)B OTFfc /(w - C) 


X (Ao(a; + gp — gp(a::)) cos kx) 


(39) 


are anti-correlated: maxima of U{x) coincide with min¬ 
ima of n(x), see Fig. 6 (d). In the strong screening case, 
the density n(x) is given, in the first approximation, by 

n{x) « d{C-U(x)) , < I. (33) 

Thus, the density profile is a mirror reflection of the ex¬ 
ternal potential, except for intervals where U{x) > C if 
they exist. In those “depletion regions” the density is al¬ 
most zero (exponentially small). To obtain the accurate 
solution for n{x), we combine Eqs. (20) and (27) into 


exp 


(^) 


+ exp 


u-C 

91T 


M - gp \ 
91T j 


= 1 , 


(34) 


which can be solved numerically^® for p at given C and 
[/, and then gp and n can be found. An example is shown 
in Fig. 3. The profile of gp(a:) [Fig. 3(a)] is seen to have 
a relatively flat minima, where screening is stronger and 
relatively sharp maxima, where screening is weaker. The 
screening is overall very efficient as the modulation of 
gp(a:) is almost an order of magnitude smaller than that 
of the external potential U{x). Since C, > Ui in this 
example, no depletion regions exist and the profile of rL{x) 
is nearly sinusoidal [Fig. 3(b)]. 

Consider now the modulation of the PL spectrum I = 
I{x,u}). Combining Eqs. (12) and (26) gives 


i{x,u}) = Bf{uj - C)Ao{uj + gp - gp(a:)). (35) 


are the Fourier amplitudes. Neglecting the broadening of 
the spectral function once again, Ap(a;) —>■ 2TrS{uj — gp). 



FIG. 6 . Illustration of the relations among spatial modula¬ 
tions of the basic quantities in the problem, (a, b) Schematic 
reproductions of Fig. 5(a, b) for one particular 7 . (c, d, e) 
Variations of the lattice potential U{x), exciton density n, 
and the exciton band edge gp over one lattice period. Note 
the reversed axes direction in (b) and (c). 
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FIG. 7. Optical resolution functions, (a) Point-spread func¬ 
tion for NA = 0.24. (b) Optical transfer function. Parameter 
Q [Eq. (36)] is the largest momentum k admitted by the opti¬ 
cal imaging system, (c) Optical transfer function as a function 
of the numerical aperture for the principal Fourier harmonic 
of momentum q = /A., A = 2 pm. 


the corresponding spectral moments are 

/r = (2-4.o)i?OTF, 

X {/{soix) - C)[eo(a;)]™cosfca;). 

Here and below {z{x)) = A“^ z{x)dx is the average of 
a given function z(x) over a lattice period. The optical 
transfer function^® OTFis the Fourier transform of the 
PSF: 


OTFfc = 


^0(Q-|k'|)0(Q-|k'-fcx|) 


2 


TT 




(41) 

(42) 


for k < 2Q. Note the normalization: OTFq = 1. At k > 
2Q, OTFfc is defined to be zero, which means that such 
harmonics are not resolved. The plot of OTFfc is shown in 
Fig. 7(b). The minimal numerical aperture necessary to 
observe at least the principal harmonic k = q is therefore 


Because of smallness of OTF^, the local PL energy is also 
approximately sinusoidal: 


E{x) 


7(1) , 7(1) 

Iq + Iq cosqx 

7(0) , r(o) 

Iq + Iq cosqx 


E{0) + Ai?avg 


SE 

cos qx. 


(47) 

Since U{x) is correlated positively with E{x) but neg¬ 
atively with I{x), the modulation amplitudes 61 and 
6E are positive. Once the electrochemical potential ( 
is known, these amplitudes can be calculated numeri¬ 
cally using fl{x) found by solving Eq. (34). We carried 
out such calculation for several sets of representative pa- 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

NA NA 




Ui (meV) Ui (meV) 


NA 


min 


g 

2ko 


(43) 


In our experiments (Sec. Ill) we are quite close to this 
lower limit, with NA « 0.24 and OTFg « 0.07, see 
Fig. 7(c). This once again confirms that accounting for 
the optical resolution effects is important. Under the 
condition 1 < NA/NAmin < 2, which is satisfied in the 
experiment, all higher harmonics of spatial modulation 
are unobservable. Hence, the spectral moments /^™^(x) 
have the sinusoidal form, e.g., 


r/ N r <5/ 

I {x) = iavg + cos qx 


I — 7^°^ 

Javg — ^0 ; 

— =40TF (/(^(O)-^(^))cosga;) 
7avg 


SI = 27f), 


{f{Ei0)-flix))) 


(44) 

(45) 

(46) 



FIG. 8 . (Color online) Parameter dependence of the PL en¬ 
ergy 5E and the normalized intensity SI/Evg, modulations. 
(a,b) SE and (57/7avg as functions of NA for T = 4.2 K, 
Ui = 4meV, 7 = 7. (c,d) SE and SI/Ia,vg as functions of 
temperature for NA = 0.4, Ui = 4meV, 7 = 7. (e,f) SE 
and SI/Evg as functions of the lattice depth for T = 4.2 K, 
NA = 0.4, 7 = 7. (g,h) SE and SI/Evg as functions of the 
correlation parameter 7 for T = 4.2 K, NA = 0.4, Ui = 4 meV. 
The average “blue shift” AEavg for each color is indicated in 
the legend of panel (b). 
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rameters, see below. In each set we fixed AE^vg, the 
directly measurable quantity. The corresponding ^ was 
determined by the standard root-searching algorithms.^® 
The results are presented in Fig. 8 . As one can see, the 
dependences of 61 and 6E on NA exhibit a threshold 
at NAmin and an approximately linear growth thereafter 
[Figs. 8 (a) and (b)]. These trends are inherited from the 
OTFq (Fig. 7). Figures 8 (e) and (f) show that SE and 
SI scale linearly with the lattice potential depth Ui when 
it is small enough. The analytical expressions in this 
linear-response regime, 


AFiavg < giT , Ui <^T + AAa 


are as follows: 

SE ~ OTFa 


UiT 


SI 


T + AAa, 


SE 


(48) 

(49) 


In this regime the modulation amplitudes are roughly 
independent of 7 ; however, they depend on the exciton 
temperature T and the ratio of measured SI /lavg and SE 
can be used to determine it. 

As Ui increases and reaches Ui ~ T-|-AAavg, the linear- 
response formulas cease to be valid. The density profile 
develops depletion regions around the maxima of U{x), 
e.g., Xq < X < K — xq where cosqx^ = 1 — {2(/Ui). 
Outside the depletion regions, for example, at |a;| < Xq, 
the density profile can be approximated by a vertically 
shifted cosine function: n{x) oc cosqx — cos ( 70 : 0 . The lo¬ 
cal energy eo{x) at |a:| < Xq is close to ( and no longer 
tracks Ui. For this reason, the energy modulation SE 
goes through a maximum and then decreases. Concomi¬ 
tantly, SI/lavg levels up at a plateau, see Figs. 8 (e) and 
(f). Under the described conditions. 


SI _ 2qxo - sin 2qxo 

lavg sin qxo — qxo cos qxo 


(50) 


which is a function of the ratio C/Uq and independent 
of T or 7 . The maximum possible SI/lavg = 4 0TFq, 
which is approximately 1.5 for NA = 0.4, is reached in 
the limit of U; —>■ 00 or T —>• 0 or 7 —0, see Figs. 8 (f), 
(d), and (h), respectively. The energy modulation SE 
vanishes in each of these three limits, see Figs. 8 (c) and 
(g). We conclude that the dependence of SI/lavg on the 
interaction parameter 7 is weak at both small and large 
Ui- However, this dependence is reasonably strong at the 
crossover point Ui ^ T+AEavg from the linear to nonlin¬ 
ear screening regime, see Fig. 8 (h). This is also the “sweet 
spot” in terms of the 7 -dependence of SE, see Fig. 8 (g). 
Similarly, it appears that the intermediate temperature 
range T ^ Ui offers the best conditions for estimating 
7 from the measured SE and SI. Below in Sec. Ill, we 
present an experimental test of this estimation method. 


III. EXPERIMENTS 

Our experiments were performed on CQW structure 
grown by molecular-beam epitaxy. Two 8 nm-wide GaAs 




FIG. 9. Sample structure and experimental results for indirect 
excitons in electrostatic lattices, (a) Schematic of the CQW 
sample (top) and its band diagram (bottom). GaAs quantum 
wells are positioned within the insulating Alo. 33 Gao. 67 As layer 
(white) surrounded by conducting n"*" GaAs layers (gray) at 
the bottom and interdigitated electrodes on the sample sur¬ 
face at the top. Labels “e” and “h” indicate an electron an a 
hole bound into an indirect exciton. (b) Electrode schematic 
(top) and the lattice potential profile (bottom) for Vb = 3V 
and SV = 0.5 V from electrostatic simulations, (c) Top: mea¬ 
sured PL intensity profile of indirect excitons. Bottom: the 
same data after subtraction of a smooth background and nor¬ 
malization, showing the periodic modulation more clearly, (d) 
Top: PL energy profile of indirect excitons. Bottom: the same 
data after subtraction of a smooth background. Experimental 
parameters for (c) and (d) are: bath temperature Tbath = 1 K, 
lattice depth Ui = 4.8 meV, photoexcitation power P = 8 pW. 


quantum wells separated by a 4nm-thick Alo. 33 Gao. 67 As 
barrier were positioned 100 nm above the n^-type GaAs 
layer within an undoped 1 /rm-thick Alo. 33 Gao. 67 As layer. 
The conducting Y GaAs layer at the bottom served as 
a ground electrode, see Fig. 9(a). Semitransparent inter¬ 
digitated top electrodes were fabricated by magnetron 
sputtering a 90 nm indium tin oxide layer. These elec¬ 
trodes generate a laterally modulated electric field per¬ 
pendicular to the quantum well plane, which couples to 
the static dipole momentum ed of the The av¬ 

erage potential was fixed by the average voltage Vb = 3 V 
and the lattice depth Ui was controlled by the voltage dif¬ 
ference SV, Fig. 9(b). The measurements were performed 
in an optical dilution refrigerator. The refrigerator had 
high stability with the vibration amplitude well below 
the lattice period A = 2.0 /im. The IXs were generated 
by a cw 633 nm laser focused to an excitation spot of di¬ 
ameter ~ 10 gm.. The exciton density was controlled by 
the laser excitation power. 

Examples of the PL energy and intensity profiles mea¬ 
sured for IXs in the lattice are presented in Fig. 9. Fig- 



























FIG. 10. (Color online) Measured and simulated modulation 
of energy and emission intensity of indirect excitons in a lat¬ 
tice. (a,b) Measured energy (a) and emission intensity (b) 
modulation for indirect excitons as a function of the lattice 
depth for different laser excitation powers (shown in the leg¬ 
end). Tbath = IK. (c,d) Simulated energy (c) and emission 
intensity (d) modulation. For each curve, the electrochemical 
potential is selected to correspond to the measured exciton 
interaction energy at the excitation power in (a,b). 7 and T 
are optimized to obtain best fit to the experimental energy 
and intensity modulation curves. 


ures 10(a) and (b) depict 6E and (5///avg as a function 
of Ui for several different excitation powers. The same 
quantities computed using the model of Sec. II are shown 
in Fig. 10(c) and (d). For each simulated curve, the elec¬ 
trochemical potential was adjusted to match the average 
PL energy in the excitation spot center. The parameter 
7 and IX temperature were optimized to obtain best fit 
to the experimental data for both energy and intensity 
modulation curves (Fig. 10a, b). 


IV. DISCUSSION 

IX gas in CQW is a model system for studying dipolar 
matter because many of its basic physical parameters can 
be controlled experimentally. For example, the range of 
accessible n can span several decades. This opens an op¬ 
portunity to experimentally test various theoretical pre¬ 
dictions regarding how many-body correlations of dipo¬ 
lar bosons evolve as a function of the particle density. 
Such a comparison with theory requires the development 
of a method for accurately determining the exciton den¬ 


sity n in absolute units that remains a challenging prob¬ 
lem. For example, estimation of n from time-integrated 
exciton emission suffers from large uncertainties.®^ The 
method for determining exciton density by measuring the 
Landau level filling factors of electrons and holes®® is ac¬ 
curate, however it requires high magnetic fields. In that 
regime the exciton properties are strongly modihed com¬ 
pared to the zero-field case. In fact, in the limit of high 
magnetic held the interaction between spatially direct ex¬ 
citons vanishes®^ that was verihed experimentally.®® The 
recently proposed technique of remote electrostatic sens¬ 
ing is promising but challenging to implement.®® Com¬ 
pared to all of the above, the lattice-based method pro¬ 
posed in our earlier work^^ and developed further in the 
present article appears to be an attractive alternative. 

The qualitative agreement of experimental and simu¬ 
lated SE (Fig. 10a, c) and SI/I (Fig. 10b, d) constitutes 
a proof-of-principle demonstration of the method. The 
accuracy of the data is still low as seen by the large data 
scattering in Fig. 11. The presented theory and experi¬ 
ment provide a guide for improvement. Since both mod¬ 
ulation amplitudes increase with NA (Fig. 8 a, b), it is 
advantageous to use an objective with a highest possible 
NA. Larger period of the lattice potential may also be 
helpful. Although the fitted values of the exciton tem¬ 
perature T (Fig. 11, inset) are in agreement with those 
from earlier studies, ®®“^®’®^^®® the accuracy of the method 
can be further improved if T is known independently. 
This can be achieved by measuring the exciton PL after 
a pulsed excitation in a time interval of the order of a 
few nanoseconds after the excitation pulse. Under such 
conditions, the IXs cool down and their temperature ap¬ 
proaches the bath temperature.^^’®® Using a defocused 
laser excitation spot to create an IX cloud with a more 
uniform density may also be beneficial. Such improve¬ 
ments of the proposed method is a subject for future 
experiments. 

On the theory side, better understanding of the disor¬ 
der effects on the exciton PL spectrum is necessary. This 
may be particularly important for correctly interpreting 
the data points for low powers [such as 1 fiW and 2 fj,W, 
black squares and red dots in Fig. 10(a) and (b)]. Here 
we confine ourselves to the following brief discussion. As 
mentioned in Sec. I, efficient loading of excitons into the 
lattice requires working with systems where the blue shift 
AUavg is higher than the characteristic energy scale of the 
disorder. In this case the interaction of excitons with ran¬ 
dom potential of disorder can be treated perturbatively. 
Thus, we can use the linear-response (effective) dielectric 
function e*(fc) to determine how the exciton gas screens 
this random potential. The zero-k limit of e*(A:) is given 
by Eqs. (30) and (31), while its large-fc behavior can be 
shown to be of the form®® e*(A:) = 1 -|- (STryn/fc^). Since 
e* depends on n, so does the disorder-induced self-energy 
correction Assuming the bare random potential is a 
white noise of bare strength Vd, we use the first Born 
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FIG. 11. (Color online) Estimated IX correlation parameter 
7 obtained from fitting the experimental data for both energy 
and intensity modulation curves in Figs. 10(a) and (b). The 
values of 7 are normalized to 7 c = 7 [Eq. (9)]. The inset 
shows the temperature values for these fits. 


approximation to get 
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r = vivd ■ 


The next improvement is the self-consistent Born approx¬ 
imation,®^ which is obtained replacing Eq in the logarithm 
with its renormalized value Eq- For T ^ T, the energies 


u! relevant for the PL are uj — Eq ^ T. The corresponding 
real part of the self-energy is 


S 


/ 

d — 


r 

7T 



(52) 


Using the definition (7) of 7 and replacing nt by AE^vg, 
we find that the random potential creates a fractional 
correction 


Jd 

7 


1 r 

TT AE^Yg 


In 



(53) 


to the correlation parameter 7 . For AE^vg = 0.5-2 meV 
in Fig. 10(a) and (b), this correction can be as high as 
~ 50% assuming r,T ~ ImeV (Figs. 4 and 11). 

Finally, we wish to point out that our method of de¬ 
termining the correlation parameter of an interacting gas 
is quite general. Similar techniques can be applied for a 
variety of other systems that can be subjected to lattice 
potentials and whose spatially modulated parameters can 
be measured. Examples include cold atomic clouds,®^’®® 
exciton polariton systems, ®^^®® and two-dimensional elec¬ 
tron layers.®^ 
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